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Abstract 

A brief review of the problem of the determination of hght-quark masses from QCD 
sum rules for the correlators of scalar and pseudoscalar currents is given. Special at- 
tention is paid to the description of the large-Nf results, obtained in collaboration with 
Broadhurst and Maxwell, for the scalar correlator, and estimates of its higher-order pertur- 
bative uncertainties within the renormalon-inspired large-/3o expansion in the MS scheme 
are given. Brief discussion is presented of the results of calculations of higher-order per- 
turbative QCD corrections to the relation between pole and MS-scheme running-quark 
masses. Their comparison with the different estimates is also given. 



1 Introduction 

The proposal to extract the values of light-quark current masses from the QCD sum rules 
for the correlators of quark scalar and pseudoscalar currents is rather old (see e.g. the works 
of Refs. 10-01 and the reviews of Refs. 0-0). However, in view of the appearance of new 
experimental data and new theoretical results and ideas, the investigations in this direction are 



continuing (for a summary of the modern situation, see forinstance the reviews of Refs. IT^, pd]]). 



When studies began, it was not clear what normalization point should be associated with the 
definite values of the light-quark current masses. Starting from the mid-90, s it became common 
wisdom, accepted by the Particle Data Group, to consider as phenomenolgical parameters the 
running quark masses of the MS scheme, normalized at either 1 GeV or 2 GeV. Another 
quantity, probably more suitable for the calculations on the lattice, is the renormalization- 
scheme-invariant mass m, connected with the running one by the following solution of the 
renormalization group equation 

r-Q^s(Q^) 7^(x) f^0o 7 7 

1^rT^^~ / 75— + 2— In 2/3c 
D /3(x) Jo /?ox Po 

= m(2/3oa.(Q2))^»/*[l + ....] , 

where the perturbative QCD /^-function and the mass-anomalous dimension function 7^ are 
now known at the 4-loop approximations (see Ref. and Refs. |T3|, Q respectively) with the 



m{Q^) = rh exp 



(1) 



one-loop scheme-independent terms /?o = [H ~ (2/3)Aj]/4 and 70 = 1. 

Another important definition of the current masses is the on-shell or pole quark masses. 
Despite the fact that, contrary to the case of charged leptons, quarks do not exist as free parti- 
cles, this definition is commonly accepted and widely used when heavy quarks are considered. 
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Within perturbation theory and the MS scheme, the running quark masses and pole quark 
masses can be straightforwardly connected. At the next-to-leading order (NLO) this relation 
was calculated in Ref. JTSf and confirmed later on in Ref. [0. Beyond the NLO level, higher 
order terms to this relation were estimated in the large- /3o approximation [|1^, |18[, which can 
be obtained from the gauge-independent large- iV/ expansion after renormalon-inspired substi- 
tution Nf — i> —Qf3o- At the NNLO level the correction to the relation between pole and running 



quark masses was estimated in Ref. within the framework of the scheme-invariant approach 



developed in the works of Ref. pO 



The results of the subsequent approximate NNLO calculations, made by evaluating expan- 
sions of the 3-loop quark self-energy for small and large external momentum and further Pade 
improvements of the corresponding mass-dependent series turned out to be in reasonable 



agreement with estimates, obtained by both large-/3o resummation [H, 18] and scheme- invariant 
methods |jl9|. In its turn, the numerical result of the approximate calculation of Ref. 
confirmed by the recent explicit analytical calculations, performed in Ref. 



JlJ] was 
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In this report the survey of the uncertainties of recent QCD sum-rule determinations of 
light-quark masses from the correlators of scalar and pseudoscalar currents at 3-loop [^]- |25] 
and 4-loop |26| perturbative QCD approximations is given. Indeed, a single-chain renormalon 
structure of the Borel images of the functions, related to the correlator of scalar quark currents 
was recently studied in collaboration with Broadhurst and Maxwell in Ref. [^. In this note 
special attention is paid to the results obtained in Ref. p8[. The higher-order perturbative QCD 



uncertainties for the studied functions are estimated with the help of the large- /3o approximation. 
The bridge between the analysis of similar questions in the case of the relation between pole 
and running quark masses, which contain the renormalon singularity |]29| , will also be 



constructed. The reason for the limitations of the scheme-invariant estimates of Ref. 
scalar channel will be clarified. 



in the 



2 Renormalon calculus 



The studies of renormalon effects (for reviews see, [0, ^) usually start from the consideration 
of the contributions of the internal single chain of fermion loops to some physical quantity D, 
say to the anomalous magnetic moment of the electron ||3^ or to the e"'"e~-annihilation Adler 
function From a perturbative point of view, these insertions form the basis of the 

large- iVj expansions, which are widely used at present. In principle, the coefficients of these 
expansions are scheme- dependent (for the discussion of the scheme dependence of renormalon 
contributions, see Ref. but in what follows this problem will not be analysed since the 



MS scheme will be essentially used. 

The next step in the renormalon calculus is the construction, from a perturbative expansion 
of the physical quantity -D(a) 



n+l 



(2) 



n>0 



2 



of the Borel transform 



and of the Borel integral 



B[D]{6) = Y.dn-T , (3) 

n>0 

D{a) = / e-^/''B[D]{6) d6 . (4) 
Jo 

The basic step in defining this integral is related to the analysis of the singularities of the 
Borel transform B[D{6)]. In the case when it has no singularities on the real positive axis 6, 
the integral of Eq. (|^) is well defined. The singularities, located on the negative axis 6, are 
called ultraviolet (UV) renormalons; after Borel resummation, they generate the sign-alternating 
perturbative QCD series with factorially increasing coefficients. The singularities, located on 
the positive axis 6 are called infrared (IR) renormalons. These are associated with the sign 
constant QCD perturbative expansions with coefficients, increasing like n\. In this case the 
Borel integral is ill-defined, but its Cauchy Principle Value (PV) can be used (for one of the 
first applications to the resummation of the renormalon- chain, perturbative contributions to 
the {g — 2)/2 of electron, in QED with coefficients growing like (+l)"'n!, see Ref. [0]). In 
fact, it is possible to show that the IR renormalon singularities can be associated to power- 
suppressed corrections of order {A/Q)^'^p\ where p is the number, related to the pole singularity, 
at the points 6 of the positive axis of the Borel plane (see e.g. Ref. []32| ). 



Up to now all discussions were made within the framework of large- iVj- expansion, which 
is well-defined in the cases of both QED and QCD. It should be stressed that, since the first 
coefficient of the QED /5-function is proportional to Nf, the large- A^/ limit of QED is identical 
to the large-/3o approximation. In the case of QCD the situation is more ambiguous. Indeed, in 
order the rewrite the large- A^j perturbative series in the form of the large-/?o approximation the 
naive nonahelianization ansatz |3^, namely the substitution Nf — >■ Nf — llNc/2=—Qj3o should 



be used. This procedure turned out to be rather useful in the analysis of both renormalon 
effects and estimates of the uncertainties due to perturbative corrections beyond the calculated 
ones (see e.g. Ref. pOj ). 



3 Definitions of the basic quantities 

Consider now the correlator of scalar quark currents Js = m{fi)ipip, namely the Green function 

UsiQ') = I I e"^%Js{x)Jsmod'x . (5) 



Its imaginary part Imll^ = 2ttsRs (1 + 0(1/ s)) is related to the spectral function Rs{s) of two 
Euclidean constructs, considered in the work of Ref. 
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The first one is defined through the dispersion relation for the first derivative of TisiQ"^) / Q"^ i 
while the second one is defined by the second derivative for IlsiQ^)- The D((5^)-function was 
used in the process of calculating the 2-loop and 3-loop QCD corrections to Rs{s) in the MS 
scheme in Refs. |^D|]-|^^. The 4-loop QCD logarithmic corrections to Us, and therefore the 
contributions to its imaginary part Rs{s), were evaluated in Ref. In the process of these 
calculations the following inhomogeneous renormalization group equation was used: 



(91og/i2 



da. 



I Tfl 



In leading order in 1/Q , and next-to- leading order in l/Nj, we have p8 

= 

dp 



TfN 



N- 



(9) 



75S 



2 + ^h{b)+0(^ 



where dp = 3, Cp = 4/3, L = ln(/i^/Q^) and b = TjNfas/iSTr). The expression of the 
coefficients Hn{L) of the large- A^/ expansion of H{L, b), namely H{L, b)=J2n>i Hn{L)b^"'^ , can 
be obtained from the all-order solution of the following relation [EM]: 



(10) 



n{n-l)H4L) = n - 4(L + 2)(7„ + 4(7„+i - 9(-irD„(L) 



E 



1 + SGd{S) 



exp( (L + -)5 



where hn are defined by the large- A^/ expansion of the /i(6)-function in the expression for 755: 

n~2 r /I 



n>l 



as 
3^ 



K + O 



N 
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and Qn are the numbers that are generated through the e-expansion of the following identity 



n>l J ^/>2 



2' - 3 



-1: 



(12) 



4 Unravelling the 5 = 1 IR renormalon knot 



The contributions of both UV and IR renormalons to II^ is given by the corresponding expres- 
sion for G{5) 



1-5 2-6 



1 ^8(l^» 



k=2 



(P- (1-5)2) 



(13) 
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It can be separated into two pieces as Gd{S) = G^{6) + G+{6). The Borel resummable UV 
renormalons are contained in the first term 



with singularities located at 5 < 0. At the large- A^j, limit it generates the sign-constant 
asjTiiptotic series. After performing a naive nonabelianization, namely after the replacement 
b — > — /3oas(/i^)/7r, one obtains sign-alternating asymptotic perturbative QCD series for Us- 

The IR renormalons are defined by G^{6), which has the following expression 

Contary to the vector case (see e.g. Ref. |^) one can observe the appearance in Eq. (|15[) of the 
first IR renormalon at 5 = 1. In accordance with the general rules of renormalon calculus, this 
IR renormalon generates the h? /Q"^ contribution to the Ds((5^)-function of Eq. (6). It should 
be stressed that this term cannot appear within the operator-product expansion machinery. 
Indeed, the first higher-twist operator, contributing to the D5(Q^)-function in the massless 
scheme comes from the gluon condensate ((as/'7r)G'^)o |^^, related to the IR renormalon pole 



at (5 = 2 (note that the perturbative QCD correction to the gluon condensate was calculated 
in Ref. Wi). 



The reason for the appearance of the 5 = 1 IR renormalon contribution to the D^-function 
is related to the fact that the dispersion relation of Eq. (6) is ill-defined. Indeed, contrary to 
the vector case, where the Ward identity allows us to fix the expression for liyiQ^) at zero 
transferred momentum, the expression for 115(0) is infinite. Therefore, the dispersion relation 
for the first derivative of IlsiQ^) / of Eq. (6) contains an ambiguity of order h? / Q^. For large 
momentum transfer this problem can be neglected. However, one should avoid the application 
of Eq. (6) in the process of long-distance studies, e.g. the determination of light-quark masses 
values. In this case it is necessary to use the dispersion relation for the second derivative of 

(see Eq. (7)), which is free from unphysical A^/Q^ term. This renormalon rediscovery of 
the advantages of the twice-differentiated Euclidean construct of Eq. (7), which was originally 



used in Ref. p| and later on in Refs. is setting its application for the extraction of 



light-quark masses on a more solid background. It should also be mentioned, that the spectral 
density Rs{s) of both Eq. (6) and Eq. (7) is also free from contributions of 5 = 1 IR renormalon. 
Moreover, since the factor 7r5/sin(7r5) of analytic continuation of the Euclidean constructs to 
the Minkowskian region is removing all single poles in Eq. (15) and transforming double poles 
in 5 into single poles, the Borel image of Rs{s) contains 5 > 2 IR renormalons only, which are 
associated with the contributions to Rs{s) of the operators Ok with dimension dk > 6. 
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5 Correlator of scalar quark currents and 
large-/3o expansion 



The large-/?o expansion in the MS scheme is widely used to generate perturbative series from 
the direct calculations of the renormalon-chain diagrams, which contain the information on 
both UV and IR renormalons structure of the related Borel transforms. Let us summarize 
the results of its apphcation for modelhng the behaviour of the perturbative expansions for 
basic quantities, related to the correlators of the scalar and pseudoscalar quark currents (see 
Ref. ll^). The quantities defined above in Eqs. (6), (7) have the following expansions: 

Rs{s) 

In leading order of the large-/?o expansion, the coefficients of these series, obtained using the 
naive nonabelianization procedure, can be presented as 

/3o""'A„+i (19) 

Po''^n+l (20) 
Po''^n+l (21) 

where the coefficients A^+i, A„_|_i and A„+i were calculated in Ref. in the MS scheme. The 
concrete results of calculations are presented in Table 1. 

Table 1: The values of coefficients in Eqs. (19), (20), (21). 



n 


An 


A„ 


A„ 


2 


7.33333 


7.33333 


5.33333 


3 


10.4696 


7.17968 


3.13622 


4 


32.6145 


8.48885 


11.6754 


5 


97.9534 


4.36402 


0.10978 


6 


503.887 


2.96849 


112.074 


7 


2194.28 


-54.2101 


-325.157 


8 


16465.8 


123.639 


3300.11 



One can see that the coefficients A„_|_i are increasing very fast. This feature is related to 
the manifestation of 5 = 1 IR renormalon contribution to Eq. (13). The suppression of this 
singularity in the expressions for Rs{s) and DgiQ"^) leads to a rather well-behaved perturbative 
series up to n = 6- loop order for Rs{s) and n = 5-loop order in the latter case, where the 



n>0 



3[m(g2)]2 1 + ^4 



n>0 



\ n>0 / 



(16) 
(17) 
(18) 



jNNA 
NNA 
-,NNA 
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corresponding coefficient A5 is remarkably small. The factorial growth of the related series 
starts to manifest itself at the level of higher-order corrections. 

Consider now the uncertainties of the large-/5o series, comparing the numerical values for the 
coefficients of Eqs. (19), (20), (21) with the explicit numbers, obtained as the result of analytical 
calculations, performed in Refs. (see Table 2). 



Table 2: Ratios of NNA estimates to exact results. 



Coeff. 


Nf = 3 


Nf = 4: 


Nf = 5 


1^2 /"2 


0.514 
0.507 
0.498 


0.496 
0.490 
0.484 


0.477 
0.472 

0.469 




0.354 
0.482 

0.764 


0.346 
0.565 

0.871 


0.339 
0.747 

1.081 



One can see that contrary to the applications of the scheme-invariant approach for estimates of 
higher-order perturbative QCD corrections in the scalar channel 
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the large-/5o approxima- 
tion gives the correct signs for the coefficients of all considered quantities. The reason for the 
failure of the scheme-invariant estimates at the 4-loop level will be clarified below. It should 
also be stressed that although the large-/5o predictions are of the same order of magnitude as 
the exactly calculable coefficients, the application of the naive nonabelianization procedure in 
the scalar channel has the tendency to underestimate the concrete perturbative terms by a 
factor of order 0.5. Keeping this in mind, and taking into account the results from Table 1, we 
arrive at the conclusion that the large-/3o approximation indicates, that the 5-loop corrections 
to the perturbative series for the Ds function and Rs{s), which enter in the determination of 
the strange-quark mass within QCD sum rules method (see Refs. 0,|2^- [^), are extremely 
small. Indeed, fixing Nf = 3 and as ~ 0.3 we obtain 



Ds = 1 + 3.67 (^^j +14.17 (^^ 
= 1 + 0.350 + 0.129 + 0.067 H 



+ 77.36 ( — 

V TT 

0.0002 



3 /Q' ^ ^ 

+ 2 X 1.26 f — 



(22) 



and 



Rs[s) = 1 + 5.667 (^j +31.864 



+ 89.156 — 
V vr 

1 + 0.541 + 0.291 + 0.078 + 0.008 . 



a. 



+ 98(^ 

TT 



(23) 



Therefore, large-/3o estimates indicate that the doubts about a serious underestimate of the 
errors of the QCD strange-quark-mass extraction from the pseudoscalar correlator at the (ct^)- 
order |]2^, formulated in Ref. using the positivity of a spectral function for a pseudoscalar 
correlator, cannot be related to the non-applicability of the perturbative QCD expansions at 
rather small energy scale. Moreover, the 0{a^) result of Ref . ^s(l GeV) = 205±19MeV is 
closer to the upper bound from Ref. namely 140 < ms(l GeV) < 254 MeV. At the energy 
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scale of 2 GeV the results of Ref. |26] are equivalent to the value Tns{2 GeV) ~ 148 MeV, 
which are higher than the current lattice results (for a review, see Ref. Note, however, 

that the current determination of light-quark masses from QCD sum rules for the scalar and 
pseudoscalar two-point functions does not touch the problem of the possible manifestation of the 
instanton contributions in the theoretical part of these sum rules (see Ref. [^). This is why the 



determination of the s-quark mass in Cabibbo-suppressed r-lepton decay characteristics, which 
seem to be free from instanton contributions, attracts more and more attention (see e.g. Ref. |4^ 
and in particular Ref. |^ ) . Note that the central value of the s-quark mass result of Ref. |^ , 
namely ms(l GeV) = 176±37oxp± 13th MeV, which is equivalent to ms(2 GeV) ~ 128 MeV, is 



lower than of the one of Ref. ||2^ and is closer to the lattice results, summarized in the review of 
Ref. |TT| . One of the reasons can be related to the approach applied in Ref. |BD[ , of resummation 
of the infinite subset of vr^-dependent effects that arize from analytical continuation of the 
mass-dependent quantities from the Euclidean to the Minkowskian region. Within the large-/3o 
approximation the special features of this procedure were studied in Ref. 



2i 



6 Pole vs running quark masses and 
the large-/3o approximation 



Another example where the naive nonabelianization procedure and the large-/5o estimates 
are working reasonably well, is the relation between pole and running quark masses. Indeed, 
this relation was recently evaluated analytically at the a^-order [^]. The results of these cal- 
culations, which confirm those obtained in Ref. |21| in a semi- analytical way, have the following 



numerical form p2| : 



M 



m{m) 



1 + 



+ 



+ 



TT 



3 V TT / V TT 

(0.6527A^? - 26.655iV/ + 190.595) 



-1.0414Ar/ + 13.4434) 



(24) 



of the QCD /3-function as p2 
M 



with as = as{m). This series can be rewritten in terms of the scheme- independent coefficients 

(25) 



m[m] 



1 + 



OL.s 



+ 



3 \ TT / V vr 

2 



(6.248^0 - 3.739) 



+ (23.497/3o^ + 6.248/5i + 1.019/3o - 29.94) 



where (3q is as defined above and (3i = [102 — (38/3) A^/]/16. It is possible now to deduce that the 
large-/?o approximation works reasonably well in this case as well, thanks to the cancellations 
between the term proportional to the jSi coefficient and the conformal /3-independent one, in 
the order correction of order to Eq. (pSj) the large-/5o approximation is working reasonably 
well [|. This feature is also illustrated in Table 3, taken from Ref. where the comparison 



wish to thank E. Gardi for discussions on questions related to this topic. 
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of the results of semi-analytic calculations of the O(a^) correction to Eq. (24) ||2T[], presented 
in the second column, are compared with the estimates of these terms, made in Ref. |]19[ using 



the scheme-invariant methods (namely the effective-charges approach (ECH) of Ref. |^T| and 
the PMS approach of Ref. |52|) and the large-/3o estimates of Ref. ITi 



Table 3: Comparison of the results of Ref. paper with estimates based on ECH, PMS 

and the large-/3o approximation for M/m{m). 



Nf 


mi 




g (ECH) 


I 


r|] (PMS) 


1 


g (large-/3o) 


2 


143(3) 


152.71 


153.76 


137.23 


3 


119(3) 


124.10 


124.89 


118.95 


4 


96(2) 


97.729 


98.259 


101.98 


5 


75(2) 


73.616 


73.903 


86.318 



One can see that the scheme-invariant methods, developed in other cases in Ref. pO[) ^^Iso give 
reasonable estimates. However, there are cases, where the comparison of the results of estimates, 
produced by scheme-invariant methods or the large /5o-approach, with the explicitly calculated 
corrections of the perturbative series under consideration, demonstrate larger ambiguities. 



7 Limitations of estimating procedures: when and why? 

It should be mentioned that unlike the case of the quahtative success of the large-/3o estimates 



in Ref. |28| of the coefficients of the perturbative series of quantities related to the Green 
function of the scalar and pseudoscalar quark currents, the application of the scheme-invariant 
estimates in this case result in the appearance of large and negative- order corrections in 
the expression for Rs{s) It is thus worth-while trying to understand the reason for the 

appearance of these suspicious numbers and, therefore, the limitations of the applications of the 



scheme- invariant procedure of Ref. PU|, used in Ref. 1113. These studies were done in Ref. 



and their main conclusions will be summarized here. 



It is worth-while recalling that the approach of Ref. pO| is constructed for the estimates of 
higher-order perturbative corrections to renormalization- group-invariant quantities of Eq. (2), 
which obey the renormalization-group equation without anomalous dimension term, with the (3- 
function defined as /5(a) = — Z]i>o A'^*'*'^ = ~/?oCt^(l + Z]j>i Qct*)) where Cj = Pi/Po- It is known 
that the method of the effective charges of Ref. prescribes to define the scheme where all 
higher-order corrections to D{a) will be nullified and the corresponding expression will have the 
following form: D{a) = doa. The new expansion parameter obeys the renormalization-group 
equation with the effective-charges /3-function, defined as 

P{~a) = -J2Pn~a'+' , (26) 

n>0 



9 



where the coefficients are scheme-invariant. The basic formula for generating the estimates 
from the scheme-invariant procedure of Ref. |^ is: 



in -I] 



Wn-^n) , (27) 



where f2„ are defined as 



Vt-i = di - ^cidi - 2d\ + 

/ 4 2 14 o 28 

1^4 = di ( C3 - -C2di + -ci(i2 + —d^ - —did2 + 4^ 

+ \{c2d2-cid^ + Mi^ . (29) 

at 5 loops. The idea that lies beyond the procedure of Ref. is that reasonable estimates of 
the higher-order coefficients dn of the perturbative series for physical quantities can be obtained 
if the coefficients /?„ of the process- dependent, effective-charges /5-functions are of the same 
sign and magnitude as the coefficients j3n of the /?-finction in the MS scheme. This assumption 
turns out to be true at the 4-loop level in QCD for the perturbative series of deep-inelastic 
scattering sum rules and the e"'"e~-annihilation D-function, in the cases Nf = 4 and Nj = 3 
in particular. Indeed the numerical analogues of the analytical expressions for the 4-loop QCD 
/3-function [|1^ have the following form 

p(Nf = 3) = -2.250a2 - 4.0000^ - 10.060a^ - 47.228a^ (30) 
p(Nf = 4) = -2.083a2 - 3.208a^ - 6.349a^ - 31.387a^ (31) 
P(Nf = 5) = -1.917a2 - 2.417a^ - 2.827a^ - 18.852a^ (32) 

with a = as/iT- It should be stressed that at Nf = 3 and Nf = 4 the effective-charges /3- 
functions of the polarized and unpolarized Bjorken sum rules and of the e"'"e~-annihilation 
D-function have a similar 4-loop behaviour, provided the explicitly unknown coefficients d^, in 
0,4 are modelled by the numerical expressions for Q^. However, this is not true in the case of 
renormalization-group-invariant quantities, related to the correlators of scalar and pseudoscalar 
currents. Indeed, for the renormalization-group-invariant quantity related to the D^-function 
of Eq. (7) as 

l^logDs(Q2^ 
2 rflogQ2 

the behaviour of the effective-charges /3-function in the 0{d^) and 0(a^)-approximations differs 
significantly from the behaviour of Eqs. (30), (31), (32) pSf. The coresponding expressions 



RniQ') = - r r:A. ' (33) 



p{Nf = 3) = -2.250a2 - 4.0000^ + 58.920a^ - 2148.5035^ (34) 
p{Nf = 4) = -2.0830^ - 3.208a^ + 53.8525^ - 1687.191a^ (35) 
(3(Nf = 5) = -1.9175^ - 2.417a3 + 49.356a^ - 1303.490a^ (36) 
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indicate the appearance of a spurious perturbative infrared fixed point. This spurious zero 
is compensated by the appearance of large and negative contributions. In the case of the 



i?£)-function, similar features were already observed at the 3-loop level in Ref. [42] and at the 



4-loop in Ref. WM both in the Minkowskian region; they were studied in detail in the Euclidean 



region in Ref. [Q. Here we conclude that these features demonstrate the limitations of the 
application of the scheme-invariant methods to the estimates of the 0(af)-corrections to the 
Ds function of Eq. (6), the Ds function of Eq. (7), and their spectral density Rs{s). 

However, the approximation of the large /3o-expansion, which we think is working reason- 
ably well in the case of the correlator of the scalar and pseudoscalar quark currents ||28|, also 



has definite limitations. The inapplicability of this method for the estimates of higher-order 



corrections to the DGLAP non-singlet kernel, studied in detail in Ref. |5^, can be explained by 
the absence of renormalon-type contributions to the anomalous dimensions of the non-singlet 
operators. A more interesting fact is related to the contradiction between the large-/9o estimates 
of the high-order perturbative corrections to the coefficient functions of = 2 and A^ = 4 non- 
singlet Mellin moments, for the F2 structure function [^, and the result of explicit calculations 
performed in Ref. ||5^ . The limitations of the predictive possibilities of the large-/3o method in 



this case are still unclear. 
Acknowledgements 

I am grateful to D.J. Broadhurst and C.J. Maxwell for the productive and pleasant col- 
laboration, which resulted in the publication of Ref. [28|. Different aspects of this work were 



discussed at the Quarks-2000 International Seminar in the talks of D.J.B. and A.L.K., and 
are described in part above. It is a pleasure to thank G. Altarelli, V.M. Braun, J. Ellis, J. 
Fischer, N.V. Krasnikov, S.A. Larin, N. Paver, A. A. Pivovarov, A. Schafer and F.J. Yndurain 
for discussions and constructive questions and comments. 

Part of the work on the problems approached in this paper was done in Russia within the 
framework of scientific projects of RFBR Grants N 99-01-00091, N 00-02-1742. 



References 

[1] A. I. Vainshtein, M. B. Voloshin, V. I. Zakharov, V. A. Novikov, L. B. Okun and M. A. Shif- 
man, Sov. J. Nucl. Phys. 27 (1978) 274. 

[2] C. Becchi, S. Narison, E. de Rafael and F. J. Yndurain, Z. Phys. C8, 335 (1981). 

[3] A. L. Kataev, N. V. Krasnikov and A. A. Pivovarov, Phys. Lett. B123 (1983) 93. 

[4] S. Narison, N. Paver, E. de Rafael and D. Treleani, Nucl. Phys. B 212 (1983) 365. 

[5] S. G. Gorishny, A. L. Kataev and S. A. Larin, Phys. Lett. B135 (1984) 457. 

[6] C. A. Dominguez and E. de Rafael, Ann. Phys. 174 (1987) 372. 



11 



7] J. Gasser and H. Leutwyler, Phys. Rep. 87 (1982) 77. 

8] A. L. Kataev, N. V. Krasnikov and A. A. Pivovarov, Preprint INR-P-0389 (1984); 
in Proc. Quarks-84 Int. Seminar, Tbilisi, 1984. 

9] S. Narison, Riv. Nuovo Cim. 10 (1987) 1. 



20 

21 
22 
23 
24 

25 
26 



S. Narison, Nucl. Phys. Proc. Suppl. 86 (2000) 242 [|hep-ph/9911454 
R. Gupta and K. Maltman, |hep-pli/010TT32[ 



T. van Ritbergen, J. A. Vermaseren and S. A. Larin, Phys. Lett. B400 (1997) 379 |Piep- 
ph/ 970 13901 



K. G. Chetyrkin, Phys. Lett. B404 (1997) 161 |[hep-ph/ 97037751 . 

J. A. Vermaseren, S. A. Larin and T. van Ritbergen, Phys. Lett. B405 (1997) 327 fhep- 



ph/9703284 



N. Gray, D. J. Broadhurst, W. Grafe and K. Schilcher, Z. Phys. C48 (1990) 673. 

J. Fleischer, F. Jegerlehner, O. V. Tarasov and O. L. Veretin, Nucl. Phys. B539 (1999) 
671 ||hep-ph/9803l93 . 



M. Beneke and V. M. Braun, Phys. Lett. B348 (1995) 513 [|hep-ph/9411229 



P. Ball, M. Beneke and V. M. Braun, Nucl. Phys. B452 (1995) 563 ||hep-ph/9502"300 



K. Philippides and A. Sirlin, Nucl. Phys. B450 (1995) 3 [|hep-ph/9503434|| . 



K. G. Chetyrkin, B. A. Kniehl and A. Sirlin, Phys. Lett. B402 (1997) 359 jhep^ 
ph/ 97032261 . 



A. L. Kataev and V. V. Starshenko, Mod. Phys. Lett. AlO (1995) 235 ||hep-ph/9502MS 



A. L. Kataev and V. V. Starshenko, Phys. Rev. D 52 (1995) 402 ||hep-ph/9412"305 
K. G. Chetyrkin and M. Steinhauser, Phys. Rev. Lett. 83 (1999) 4001 ||hep-ph/9907509 



K. Melnikov and T. v. Ritbergen, Phys. Lett. B482 (2000) 99 ||hep-ph/9912391|] . 

J. Bijnens, J. Prades and E. de Rafael, Phys. Lett. B348 (1995) 226 [|hep-ph/9411285|] . 



K. G. Chetyrkin, C. A. Dominguez, D. Pirjol and K. Schilcher, Phys. Rev. D 51 (1995) 
5090 ||hep-ph/ 94093711] . 



M. Jamin and M. Munz, Z. Phys. C66 (1995) 633 ||hep-ph/9409"335|] . 

K. G. Chetyrkin, D. Pirjol and K. Schilcher, Phys. Lett. B404 (1997) 337 jE^ 



ph/96l2394| 



[27] P. Colangelo, F. De Fazio, G. NarduUi and N. Paver, Phys. Lett. B 408 (1997) 340 Eep 



ph/ 9704249 



12 



[28 



[29 

[30 
[31 

[32 

[33 
[34; 
[35 
[36 
[37 
[38 
[39 
[40 
[41 

[42 

[43 
[44; 
[45 
[46 

[47; 

[48 
[49 



D. J. Broadhurst, A. L. Kataev and C. J. Maxwell, Nucl. Phys. B592 (2001) 247 jE^ 
ph/0007152| 



I. I. Bigi, M. A. Shifman, N. G. Uraltsev and A. I. Vainshtein, Phys. Rev. D 50 (1994) 
2234 ||hep-pli/ 94023511 . 



M. Beneke and V. M. Braun, Nucl. Phys. B426 (1994) 301 ||hep-ph/9402"364 |. 



G. Altarelli, CERN-TH/95-309; in Proceedings of Int. School of Subnuclear Physics, July 
1995, Erice, Italy. 



M. Beneke, Phys. Rep. 317 (1999) 1 ||hep-ph/9807"443 



M. Beneke and V. M. Braun, |hep-ph/ 00 10208 . 

B. Lautrup, Phys. Lett. B69 (1977) 109. 

V. I. Zakharov, Nucl. Phys. B385 (1992) 452. 

L. S. Brown, L. G. Yaffe and C. Zhai, Phys. Rev. D 46 (1992) 4712 ||hep-ph/9205"2Tl . 
N. V. Krasnikov and A. A. Pivovarov, Mod. Phys. Lett. All (1996) 835 ||hep-ph/960227l 
A. S. Ilchev and V. K. Mitrjushkin, J. Phys. G G7 (1981) L221. 
D. J. Broadhurst and A. G. Grozin, Phys. Rev. D 52 (1995) 4082 ||hep-ph/9410"240|] . 

C. N. Lovett-Turner and C. J. Maxwell, Nucl. Phys. B452 (1995) 188 ||hep-ph/9505^ . 



S. G. Gorishny, A. L. Kataev and S. A. Larin, Sov. J. Nucl. Phys. 40 (1984) 329. 

S. G. Gorishny, A. L. Kataev, S. A. Larin and L. R. Surguladze, Mod. Phys. Lett. A5 
(1990) 2703. 

S. G. Gorishny, A. L. Kataev, S. A. Larin and L. R. Surguladze, Phys. Rev. D 43 (1991) 
1633. 



K. G. Chetyrkin, Phys. Lett. B390 (1997) 309 ||hep-ph/96083T8[ . 

D. J. Broadhurst, Z. Phys. C58 (1993) 339. 

M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B147 (1979) 385. 
L. R. Surguladze and F. V. Tkachov, Nucl. Phys. B331 (1990) 35. 
F. J. Yndurain, Nucl. Phys. B517 (1998) 324 [|hep-ph/9708300| . 

E. Gabrielli and P. Nason, Phys. Lett. B313 (1993) 430. 

K. G. Chetyrkin, J. H. Kuhn and A. A. Pivovarov, Nucl. Phys. B533 (1998) 473 jE^ 



ph/9805335| . 



13 



[50] J. G. Korner, F. Krajewski and A. A. Pivovarov, hep-ph/0003165 . 
[51] G. Grunberg, Phys. Rev. D 29 (1984) 2315. 
[52] P. M. Stevenson, Phys. Rev. D 23 (1981) 2916. 

[53] E. Gardi and M. Karliner, Nucl. Phys. B 529 (1998) 383 ||hep-ph/980m^ 
[54] S. V. Mikhailov, Phys. Lett. B 431 (1998) 387 ||hep-ph/9804"26^ . 
[55] L. Mankiewicz, M. Maul and E. Stein, Phys. Lett. B 404 (1997) 345 ||hep-ph/9703"35^] 
[56] S. A. Larin, T. van Ritbergen and J. A. Vermaseren, Nucl. Phys. B 427 (1994) 41. 



14 



